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ABSTRACT 


The purpose of this study was the development of a computer system 
for research and instruction in signal processing for underwater sound 
applications. 

The nature of the problem demanded the use of three computers inter- 
faced by hardware and software to operate as a single system. A general 
purpose analog computer was employed to generate signals and noise, 
provide analog processing and filtering and to allow the input of actual 
Signals collected in the field. A general purpose digital computer 
provided for the examination of existing processing methods and for the 
investigation of new techniques. A high-level graphics computer was 
used for data display and as the system control station which provided 
for a high degree of human operator processing and interaction. 

All the hardware required was provided in the Computer Laboratory 
of the Naval Postgraduate School, Department of Electrical Engineering 
and this study prepared the software for its convenient and efficient 


use for this application. 
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I, INTRODUCTION 


A. PURPOSE 

The purpose of this work was to develop a convenient method and 
efficient programs to employ the unique computer system of the Naval 
Postgraduate School Computer Laboratory to build a hybrid system for 
Signal analysis and processing. The primary good of the system was to 
achieve a high speed computational tool with a high degree of human 
operator interaction. In addition the system was to provide the flexi- 
bility and interaction which would permit the operator to-use his 
heuristic capabilities in approximating solutions in a design procedure 
when the search for adequate patterns and processors is the goal, or to 
permit the student to verify theoretical results for a wide variety of 
Situations that he may configure almost at will. 

It was also assumed from the beginning that the system should be 
directly aimed to the underwater sound detection problem. This decision 
not only dictated the name of the system namely a Hybrid Sonar Simulation, 
but also directed the choice of the preprogramed applications initially 
included. This choice does not, however, restrict the application of the 
program to a much wider field of signal processing limited primarily only 


by the user's imagination and expertise. 


B. SYSTEM REQUIREMENTS 

After establishing the general purpose, some specifications were set 
up to further specify the scope to be achieved. 

In order to enable the user an easy and rapid perception of the 


Signals and the results of their analysis or processing, a graphical 








display was chosen. To allow him to interact with the process, it was 
necessary to make available means of controlling the signals, processors 
and operations at the same location where the situation was being dis- 
played and while it was shown. The fundamental requirement for this type 
Of operation is very high speed computation. 

Since comparison of signals and patterns was one of the main objec- 
tives, the user should be able to choose among different classes of 
Signals and then vary their parameters. Also the characteristics of 
filters should be under control. 

For analysis and processing there was established the requirement 
for spectral analysis, auto- and cross-correlation and convolution. 

Three types of processors, a linear cross-correlator and two non- 
linear processors were included as well as the capabilities to determine 
Signal-to-noise ratio, input-output performance characteristics ("proces- 
sing gain") and the receiver operating characteristics (ROC). 

Random noise was necessary to be mixed with the uncorrupted signals 
in controlled amounts. A Gaussian-like noise was considered adequate 
for initial test purposes. 

Finally a possibility for taking hard copies of the curves displayed 
had to be incorporated to allow record keeping for further studies and 


comparisons. 


C. ELEMENTS AVAILABLE 

The elements available in the computer laboratory and used to compose 
the system include (see Fig. 1.1): 

Digital computer — XDS 9300 — a medium sized general purpose computer 
with 32K words of main memory and a processor capable of handling a large 


set of operations with very good efficiency. In the installation it is 
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connected to the analog computer and the graphics digital computer which 
allows the hybrid configuration. Also available are several of the 
Standard input-output devices among them a magnetic drum secondary 
Storage with a two million character capacity. 

Analog computer — The Comcor CI 5000 is a very well equipped analog 
computer which includes a parallel logic portion that allows extensive 
decision making and control functions to be performed. 

Graphics computer — The AGT-10 is a small general purpose computer, 
with 8K words main memory which has as a special feature a CRI output and 
the necessary hardware to write on the screen, generate vectors and 
refresh the pictures. 

Converters — Both a high precision analog to digital converter, or 
Sampler, and digital to analog converters are available. 

Software — The system software includes the normal mathematical func- 
tions and subroutines to perform some FORTRAN operations, a set of programs 
to control the analog computer from the digital, programs to control the 
graphics terminal, two subroutines to plot curves with the line printer 
and a program to compute fast Fourier transforms according to the so called 


Cooley-Tuckey method [Ref. 2]. 


D. GENERAL APPROACH TO THE SOLUTION 
The logical approach was to use the elements available according to 
their best performance characteristics as applied to the problem. 
Naturally, the graphics computer was taken as the data presentation 
component and then became also the system control unit. Sitting in front 
of it the engineer incorporates himself into the system as the center of 


decisions. He selects an operation, analyzes the results, determines 








the next step, compares the new results, and introduces new parameters 
until satisfied with the information accumulated. Since the graphics 
terminal is itself a general purpose computer, the types of data to be 
displayed as well as the format of the display may be varied to investi- 
gate the effect on Human Operator performance. 

In nature the great majority of data or information comes in a 
continuous or analog form, the discretization being an human abstraction 
or simplification for easier handling. In particular, that is the situ- 
ation encountered in the detection problem where the signals are 
continuous. For these reasons the analog computer was utilized to 
generate the signals. The analog computer offers the possibility of 
wider range of variations with very simple and easy to perform change of 
parameters. By using signals of the same nature as the real ones, 
practically no changes are required if instead of simulated functions 
true signals are the object of analysis either in the form of a pre- 
recorded reproduction or as the output of a receiver. 

To analyze or process the signals the discretization becomes a 
necessity since the digital machines are by far more precise and very 
well adapted to the computations described by fixed algorithms. Hence 
the choice of the digital computer for this part of the system. 

Two ways may be used to show the arrangement, From a functional stand- 
point the operations and their inter-relationship are the constituent 
blocks as presented in Fig. 1.2. Depending on the manner these functions 
are employed, two modes of operation result and may be named ANALYSIS and 
SIMULATION. In the ANALYSIS mode it was intended that basically one func- 
tion be performed at a time under operator control thereby giving him an 


opportunity to verify intermediate results and to proceed accordingly. 
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Through the SIMULATION mode some pre-programmed problems were configured 
and automatic sequences could be entered to determine processing gain, 
receiver operating characteristics and detection distance. Some of the 
processing functions could also be executed by a program that performs 
correlation by a point by point multiplication and addition; since these 
Operations serve no other purpose than to demonstrate the long time re- 
quired by this direct approach they were thought of as belonging to a 
DEMONSTRATION mode. 

From a physical point of view the system is described as depicted in 
Fig. 1.3. The digital computer, besides providing most of the processing 
functions, exerted effective control of the other units from which it 


received interrupts to initiate or synchronize operations. 


E, ELEMENTS CREATED 

To implement this system several elements were created to execute the 
functions and satisfy the requirements established. A brief description 
of these elements is presented in the following paragraphs. 

ia control 

A program was written to read the control codes typed on the 
graphics terminal and compare them against a table to establish pointers 
that indicated the proper sequence to be executed. 

Another program was introduced to allow the operator full control 
of the analog computer from the graphics terminal. Through this feature 
it is possible to select modes, set or test lines, set potentiometers or 
scan the various analog components by typing simple messages at the 
console, 

2. Signal Generation 
Signal generation was developed from sine generators patched on 


the analog board. Bearing in mind the sonar application, one of these 
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generators was associated with additional analog amd logic circuitry to 
permit simulation of CW, FM-slide or FSK codes. Timis generator was also 
equipped with positive feedback to make possible a more realistic effec- 
tive sampling rate which required the analog computer to run at a very 
high frequency. Without some regeneration in this case the wave damping 
would be intolerable. 

Three other sine generators were implemented to simulate propeller 
Signatures, which present a low-pass line spectra composed of a funda- 
mental and several of its harmonics [Ref. 1]. To achieve this result the 
Sine waves were fed to a combination of multipliers whose outputs were 
Summed in controlled amounts to give a final pattern with up to six 
harmonics, each of which with any desired intensity. 

3. Noise Generation 

To simulate a random-like noise a pseudo-random binary sequence 
was generated by a maximal length feedback-shift register and then trans- 
formed into a pseudo-random analog square wave. This signal, after proper 
filtering presents Gaussian statistics at least up to a first-order [Ref. 6]. 

4. Acquisition 

Signal acquisition to the processor was normally made after filter- 
ing and then through the sampler or Analog-to-Digital converter (ADK). 

Two filters were included: an RC lowpass and a narrowband, high-Q, 
second-order filter. The latter offered the advantage of being easily 
controlled from the graphics terminal to adjust frequency, bandwidth or 
gain. 

The purpose of the filters was to represent transducer and part 
of preprocessor characteristics, but they could also simulate or be pre-~ 
ceded by other filters with any other objective, such as transmission 


medium characteristics. 
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The sampler was controlled by a program that not only synchronized 
the operation but also simulated the remaining parts of a pre-processor 
for the non-linear cases, and stored the samples, conveniently normalized, 
into the proper locations of the processor memory. 

The samples were distinguished and stored as belonging either 
to the REPLICA or to the SIGNAL. While the choice of names relects the 
application to sonar, for a general case these names serve only to identify 
independent storage areas where signals could be stored for further opera- 
tion or retrieval. Similarly, a third area where the results of correlations 
were placed, received the name of CORR. 

5. Spectral Analysis 

A special subroutine was prepared to make use of the Fast Fourier 

transform program and prepare the results for display. 
6. Processing 

A replica correlator simulation was employed. Depending on how 
the stgnals have been preprocessed the processor became a linear, 
amplitude-limited or a clipper-limited type. 

The program written to execute the correlation was based on the 
properties of discrete Fourier transforms and on the time domain convolu- 
tion theorem as applied to the discrete case. For cross-correlating the 
replica against a very long time series, as is the case where processing 
the incoming signal in a detection problem, a partition technique was 
adopted [Ref. 4]. 

The correlation could also be evaluated by multiplying both 
functions, adding the results, shifting and repeating these operations 
until the replica had been slid entirely across the signal. The time 
required for this procedure was incompatible with the real-time specifi- 


cation and the method was left for illustration only. 
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7. Performance Determination 

Under this general classification were considered the input-output 
Signal-to-noise ratio characteristic, the receiver operating characteris- 
tic and distance measurement. 

A program was prepared to automatically sample the replica and 
the signal, form the signal either with noise alone or with noise plus 
transmission modulation code, process the received samples, compute power, 
determine threshold and threshold crossing, evaluate probability of 
detection or false-alarm and measure distance, all in the proper sequence 
according to the case. 

Power determination, when necessary, was computed through the 
application of the Wiener-Kinchin relation and Parseval‘s theorem. 

The threshold was normalized to the peak value of the correlator 
Output when no noise was present. 

When using these automatic operations, where only the final 
results and not the waveforms were of interest, the sampling rate utilized 
was related to twice the information -Eandwidth. This allowed processing 
samples with a more realistic information content. Attempts to display 
the waveforms at these sampling rates showed that while the signals were 
completely distorted, they revealed amplitude variations at a rate 
comparable to the significant bandpass. For signal analysis, when the 
waveforms were required to be displayed without great distortion, the 
sampling rate had to be considerably greater than twice the highest 
Frequency present. 

8. Display 
Programs were prepared to display the functions and results 


covering all cases treated. 
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Figure 1.4. Spectrum of a CW signal. This figure is a photograph of 

the actual display. It shows the time function, absolute value and 

phase angle of the transform. The location of the impulses, at 15000 Hz, 
serve as reference for the frequency plots in the next pictures. 
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Figure 1.5. Spectrum of an FM signal formed by a 5000 Hz sinu- 
soidal linearly modulated by 1000 Hz. Comparison with Fig. 1.4 
Shows the wider band covered. Discrepancies observed in the phase 
angle plot here, as in some of the following pictures, are due to 
digital computer noise in the computation of arctan. 
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Figure 1.6. epee of a sinusoidal signal of 5000 Hz two-tone 
frequency modulated by £500 Hz. 
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Figure 1.7. Spectrum of unfiltered noise. 
1s a pseudo-random square wave wh 
a low-pass filter with cut-off frequency beyond the region of 
materest. This filtering is necessary to assure Gaussian first- 
Order statistics to the simulated noise. In the band up to 5000 Hz 
(compare with Fig. 1.4) no significant peaks appear in the spectrum. 


The unfiltered noise 
ich has been passed only through 
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9, Copyin 
A subroutine was written to copy the data displayed on the screen 
into the format required to use the two programs for line printer plotting. 
Also a sequence in this subroutine was prepared to send the 
information to the analog computer, through the digital to analog 
converter, in order to plot the curves with a paper recorder, This tech- 


nique proved very useful in several applications. 


F. SYSTEM CAPABILITIES 

The system was made flexible enough to allow a wide range of appli- 
cations. Some of these were checked with very impressive results as 
demonstrated below. 

Initially, Figures 1.4 through 1.8 show types of signals used. 
Fig. 1.4 presents a CW burst of 5000 Hz and 2.53 msec; also shown is its 
Spectrum, magnitude and phase. In Fig. 1.5 the 5 kHz signal is linearly 
FM modulated by 1000 Hz and in Fig. 1.6 the 5 kHz carrier has an FSK 
modulation of £500 Hz. Fig. 1.7 depicts a 2.53 msec sample of unfiltered 
noise; the frequency plot covers the range of t50 kHz. Fig. 1.8 is a 
simulated propeller signature with a fundamental of 60 Hz and six harmonics, 

Figure 1.9 shows an example of filter characteristics determination 
through the application of an impulse at its input. In the upper part the 
impulse response is presented followed by the magnitude and phase of the 
transfer function. The filter tested was a narrowband centered at 5 kHz 
with: a bandwidth of less than 1 kHz. 

In the next figure, Fig. 1.10, the same random noise previously dis- 
played in Fig. 1.7 is seen but after being passed through the narrow 
bandpass filter. This experiment was a demonstration that noise behaves 


as a sinusoidal signal when narrow-band filtered. 
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Figure 1.8. Representation of a low<pass periodic signal Similar 

jin characteristics to a propeller-generated noise. The line spectrum 
shows well defined harmonics up to the sixth. The fundamental 
frequency was 60 Hz. 
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Figure 1.9. This figure presents the narrow-band filter impulse 
response, and both the amplitude and phase angle of its transform 
or transfer function. 
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Figure 1.10. In this picture the noise spectrum is seen again but 
now after narrow-band filtering. The spectral density is similar 
to that of a sinusoidal signal but the phase angle has a random 


variation. 
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Figure 1.11. The FM-signal of Fig. 1.5 is shown after passed 
through the narrow-band filter. The effect of the filter 
narrow-pass is reflected in the time function by a reduction 

in amplitude for the frequencies above and below the center 
frequency. In the frequency domain a reduction in bandwidth is 
observed. The sample in this example is twice as large as that 
of Fig. 1.5 and shows the effect of the transient filter response. 
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Figure 1.12. This picture shows how the filter output of the 
previous figure can be obtained by computation. Above is displayed 
the filter impulse response time-reversed. In the middle the 


unfiltered FM-signal, and below the result of the convolution of 
these two functions. | 
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Figure 1.13. Power spectral density of the signal presented in 
Fig. 1.8 but now highly contaminated by noise, The auto-correlation 
function and its transform are shown. 
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Figure 1.14. Auto-correlation function of a CW<pulse. 
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Figure 1.15. Auto-correlation of the FM signal shown previously 
in Fig. 1.5. Comparison with the auto-correlation of the CW 
Signal presented in Fig, 1.15 shows the better definition of the 
correlation as a consequence of the wider bandwidth. 
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Figure 1.16. Auto-correlation of the FSK pulse presented in Fig. 1.6. 
By comparison with Fig. 1.15 it is seen that the uneven distribution 
of the energy over the bandwidth affects the result of the correlation. 
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Figure 1.17. Auto<correlation of unfiltered noise. The picture 
clearly shows the lack of correlation among random-noise samples. 
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Figure 1.18. In this picture, the CW signal of Fig. 1,4 was 
amplified and clipped after having been passed through a narrow- 
band filter. The envelope of the auto-correlation is similar to 
that of the undistorted signal, which shows that narrow-band 
Signals may be treated by absolute-value processor models. 
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In Fig. 1.11 the FM signal of Fig. 1.5 is again displayed but now 
after filtering. In Fig. 1.12 the same result is obtained by convolving 
the Fier impulse response with the unfiltered signal — the upper part 
of the figure is the filter response flipped, the central the unfiltered 
Signal and the lower the result of convolution. This experiment demon- 
Strated the capability of the system to simulate digital filters or any 
other filter once the impulse response is known. 

Figure 1.13 shows the result of power spectral density evaluation of 
the propeller signature already mentioned but now noise added and after 
lowpass filtering. Displayed are the auto-correlation of the signal and 
the magnitude and phase of the power spectrum. 

The next four figures, 1.14, 1.15, 1.16 and 1.17 depict the auto- 
correlation of the signals previously shown in Figures 1.4 through 1.7, 
i.e€., CW, FM, FSK and noise. 

Figure 1.18 shows the auto-correlation of a CW signal processed by a 
clipper-limited type of processor. Since only the zero crossing informa- 
tion is kept, the resultant signal is a binary signal. 

The following figures present the results of the cross-correlation 
of a replica and a signal composed of the transmission pattern plus noise. 
The cross-correlation was normalized in such a way that if the transmis- 
Sion code only was present the peak output would be one. All samples are 
taken after narrowband filtering and correspond to the FM-slide signal. 

Figures 1.19 and 1.20 refer to a linear cross-correlator output for 
two different signal-to-noise ratios while Figs. 1.2] and 1.22 use the 
Same SNR but applied to a clipper-limited processor. The experiment 
showed that for high signal-to-noise ratio the non-linear processor has 
a much lower sensitivity than the linear one; the differance, however, 


is almost non-existent for a low-input signal-to-noise ratio. 
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Figure 1.19. Cross-correlation of an FM-replica against this 
signal plus noise at a high signal-to-noise ratio. The periodic 
or circular nature of the Discrete Fourier Transform is evident 


from the result of the correlation. 


34 


il al 


“ANS 


fuse]! 


tee \ Aah Rua 
i a \s al. 


Y i ihe 
i he 


(1 | : 
pO N EN gy ile o~ wll 


\ 
hi Te si “1 i A pants 7 
eens ea rank WY fie 
( Wirral why! Wile 
a ti, 
\ 


TORR CROSSCORRELATIO“N REPLICA: PM-GUIDE SIDI: PM-GiidE AUS 'OISE GR iB 
FERICS vee ia ENTER 1 LStDER FGRST LETTEP 





Figure 1.20. This picture is the result of the same operation as 
that used to obtain Fig. 1.99. but applied when the signal-to-noise 
ratio is low. Although the functions are normalized to the same 
maximum amplitude in the display, the deterioration caused by 
noise is evident, especially when the central "bulge" is observed. 
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Figure 1.21. This figure shows the result of processing the 
Signals of Fig. 1.19 with a clipper-limited processor. 
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Figure 1,22. This figure presents the results of applying a 
clipper-limited processor to the signals of Fig. 1.20, where a 
low signal-to-noise ratio was employed. Since the amplitudes 
are displayed normalized to the same value, comparison can be 
made by considering the peak-to-peak oscillations of the outer 
segments. At low SNR the non-linear processor approaches the 
linear processor better than when the SNR jis high. 


3) 








For the next experiment (Fig. 1.23) in order to simulate a doppler- 
shifted reception, the FM signal was modified, but not the replica. The 
input signal-to-noise ratio is the same as for the case of Fig. 1.19 in 
order to present the degradation caused by doppler. 

In Fig. 1.24 a replica of a propeller Signature is cross-correlated 
against a long sequence of incoming samples where the desired signal is 
present but at a very low signal-to-noise ratio. The replica was corre- 
lated against the long sequence using the partition method; the result 
displayed is only a section of the correlation corresponding to 9.12 
msec. It was observed that the peak output was low, although the oscil- 
latory character of the correlation seemed to carry much of the information 
required for detection. A Fourier transform of the result, illustrated 
in Fig. 1.25, showed how this information could be identified. 

The last figure of this series, Fig. 1.26 is the result of an experi- 
ment similar to that of Fig. 1.20 but with the sampling rate just greater 
than twice the bandpass (up to this point the sampling rate used was 
twenty times the highest frequency). It was noted that the displayed 
replica is completely distorted (no narrow band filters were used to 
recover the signal from the samples) but revealing an amplitude variation 
in the order of the bandwidth. The cross-correlation, although of unrecog- 
nizable shape, has a highest peak. In this case, where both the total 
number of samples and the bandwidth have been conserved, the reduction in 
the sampling rate corresponded to a larger integration time with less 
redundant samples being processed. Conversely, in the case of Fig. 1.20 
the integration time was short and several non-independent samples were 
processed, contributing nothing to the result. 

The importance of input-output performance and ROC determination for 
comparison of pracessors and transmission codes was verified by several 


problems simulated as described in a later section, 
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Figure 1.23. This figure presents the result of an experiment 
designed to simulate a doppler shift of the incoming signal. 
Comparison with Fig. 1.19 demonstrates the degradation in the 
correlation as a consequence of the shift. 
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Figure 1.24. In this picture a section of the cross-correlation 
of the propeller noise replica and the continuously received 
Signal with large amount of noise is shown. Note that the 
cross-correlation function, if not further processed, does not 
reveal any information. 
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Figure 1.25. This figure shows the result of the spectral analysis 

of the cross-correlation function displayed in Fig. 1.24 

Comparison with Fig. 1.8 demonstrates that the presence of the 

original signal is now readily identifiable. The same is not true, 
eet 7 the power spectrum of the signal before correlation 
mig. ais) 
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Figure 1.26. This figure presents the effect of sampling the 
Signals of Fig. 1.20 at a sampling rate close to twice the signal 
bandwidth. Note that while the signal is badly distorted at this 
low sampling rate, and without any reconstruction filtering, the 
peak of the cross-correlation is much better defined. 
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G. DETAILED DESCRIPTION 

In the following sections a detailed description of the system imple- 
mentation together with a review of the mathematical concepts utilized 
will be given. Further, a brief section on the graphics terminal and 
discussion of sample simulations is presented. Finally, some conclusions 
and recommendations for additional applications are provided, 

The author does not claim originality for the theory discussed, which 
can be found in good reference books such as those mentioned in the text. 
Perhaps only a certain degree of innovation exists in the interpretation 
given to information bandwidth which seems to depart somehow from the 
conventional. 

The work described here was, and was intended to be, an engineering 
study. A program requiring a solid theoretical background, a clear view 
of the objectives, and a willingness to make the unavoidable compromises 
characteristic of real problems, in order to transform the resources 


available into a useful product. 
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ITI, DIGITAL PROCESSING 


The system was implemented to process the signals digitally, since 
the great majority of processors of interest are digital. The reason for 
that rests mainly upon the ability of digital machines to handle algorith- 
mic operations very rapidly and with high precision. The discretization 
of data required for digital processing does not necessarily represent a 
reduction in information. If the sample rate is properly adjusted, only 
independent samples may be processed and no time is wasted with the 
redundancy existing in the original signal. 

To use the XDS-9300 of the laboratory it was necessary to write 
programs in order to simulate the system and processors. Since this is 
a general purpose computer several other processors can be simulated 
without any hardware change, requiring only the addition of new programs 
or modification of those initially included. 

In this section a description of the computer program, its usage, and 
the theoretical reasons determining its construction are presented together, 


in order to make it easier for the reader to follow the discussion. 


A. OVERLAY STRUCTURE 

Since, as noted before, the computer has to store in core a large 
collection of library subroutines to be used in the program, including 
those to control the analog computer and the graphics terminal, and also 
has to store relatively large arrays for the computations of correlations 
and convolutions, little space is left to the program itself. To overcome 


this difficulty the program was segmented as shown in Fig. 2.1. 
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In brief, the purpose of each segment is as follows. Segment MAIN 
contains the program initialization and the pivot points where the status 
is saved when a new overlay is called in. Segment G is composed of the 
subroutines required to control the analog computer from the graphics 
terminal. Segment A permits bringing in functions and performing opera- 
tions on them under operator control, and constitutes the base of the 
ANALYSIS mode. Segment AB is the principal part of the processor and 
computes the cross-correlation of two functions; it may be called either 
from ANALYSIS or PROCESS through a bridge in the main segment. Segment 
P allows the representation of a processor to be automatically put to 
work and enters what is called the SIMULATION mode. Segment T is designed 
to copy what is presented on the screen. Segment D computes correlations 
through a direct approach for a DEMONSTRATION mode. Finally, segment MA 
provides graph and test initialization, decodes the instructions typed 
in by the operator, and selects the path to be followed. This segment 
provides the capability to select the type of processor or pre-defined 
Signals, and to enter new parameter values. 

If the system is to be modified by addition of new segments, the 
maximum size allowed is approximately 4500 (octal) words per segment. 
This limit however depends on the addition or deletion of library sub- 


routines and also on the size of the sample as will be noted later. 


B. MAIN PROGRAM 

The function of the main program, besides serving as a bridge between 
the various operations and providing initialization of the simulation, is 
to define permanent Tocation in memory for some constants and the arrays 


where the samples and results of operations are stored. 
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The constants are stored in the COMMON area and their meaning will 
be explained with the subroutines in the Appendix. One exception to this 
rule is the constant M that defines the sample size; it can be modified 
by the operator and is initially set to 256, its maximum allowable value. 
The value of M must be a power of 2 due to the methods used for computation 
of correlation and Fourier transforms. 

The arrays defined in the main program are passed to the subroutines 
as arguments; with this approach it is readily possible to make use of 
the same storage area for different applications according to the situ- 
ation, and it is not necessary to modify the subroutines when a new value 
is used. 

The sample size, M, serves as a parameter to define the blocks of 
array storage. At this point it is convenient to remember that in the 
XDS-9300 a real variable uses two words of storage while an integer 
requires only one. It follows then, that one block of M locations actually 
requires 2M words. 

The map of Fig. 2.2 shows the correspondence of designations through- 
out the program. 

The basic assignment of the arrays is as follows: 

YR — replica, or the signal transmitted; size M 


YS — signal, or received signal; size 3M, what allows a faster 
computation of recursive correlation 


YC — correlation or result; size 2M 
VAGO — empty buffer used for several applications; size 5M 


IB — buffer initially designed for text and graph blocks; size 2M 
(plus one word) 


Before any sample is entered, the main program creates a dummy replica 


and signal, two different size rectangular pulses which allow easy 


verification of the system functioning. These two functions are immediately 
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FIGURE 2.2 
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dumped into secondary storage. This is also always done with the replica, 
Signal and correlation in the ANALYSIS mode. However, only the very last 


one of each of these functions entered into the system is saved. 


C. SUBROUTINES SEGSAMP, XEGSAMP AND SAMPLE 

Subroutines SEGSAMP and XEGSAMP have the responsibility of controlling 
the sampling of the analog computer. They are identical but have distinct 
identifiers since they appear in two different levels of segmentation. 
They also provide normalization so that the amplitude of the samples is 
less than a maximum value compatible with the system. 

Another function of SEGSAMP is to prepare the signals for the simula- 
tion of three types of processors: linear cross-correlator, amplitude- 
limited cross-correlator and clipper-limited cross-correlator. 

SEGSAMP controls three types of sampling according to the value of 
INDEX, 1.e., replica, initial signal, or continuation signal. 

The replica normalization factor, RCOF, signal normalization factor, 
SCOF, and their ratio Q are saved in common storage, because they are 
necessary to the computation of the input signal-to-noise ratio, SNR; 

Q is also necessary to correct the computed threshold value for detection. 

The processor type is determined according to ISAMP, i.e., linear, 
Clipper-limited or amplitude-1limited. 

INDEX has tne initial value of one and ISAMP is zero; these numbers 
however may be set as desired by the operator. 

The actual samples are not taken by SEGSAMP. After the subroutine 
is entered and the proper sampling sequence is selected, the computer 
is trapped in a loop until an interrupt is generated by the sample clock. 


At this moment the subroutine SAMPLE is connected and samples the analog 
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computer, storing the value into ARRAY; at the same time it decreases a 
counter, MAC, which has the effect of opening the loop in SEGSAMP and 
allows copying ARRAY into the proper location. 

The sample is taken at the circuit node connected to the trunk line 
determined by LINE whose inittial value is one, but can also be adjusted 
to any legal trunk line designation. 

Figure 2.3 is a simplified flow chart that illustrates the operation 


of the subroutines. 


D. SUBROUTINES DISPLAY AND TISPLAY 

These are also the same subroutine with different calling names to 
allow their location in two levels of segmentation. 

DISPLAY is designed to encode the functions as required for display 
on the graphics terminal. It has basically three options. The first 
option was developed for the presentation of correiation and convolution 
operations. In this case the y-axis is positioned at the left of the 
screen. The signals to be correlated and the resuit of the operation are 
plotted in three sets of coordinate axes one below the other. As a 
variation of this option only the first function and the three coordinate 
axes or only the other two functions are displayed at one time; this 
feature is useful to the subroutine FASTCOR as is shown later. The number 
of points used to describe each function may also be varied. The second 
option is applied to the Fourier transform and here the same number of 
points are required for all three functions. The y-axis is placed at the 
center of the screen. In this case the first one is the time function 
itself, the second the absolute value of the transform, and the third the 
phase angle. The last option is pravided for single function display. 
The function is presented at the center of the screen with the ordinate 


expanded, 
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E. SUBROUTINE FASTCOQR 

This subroutine constitutes the basic feature of the system and its 
epose is to perform a fast convolution. Next will be presented a review 
of the mathematical concepts that were considered to build this program. 


The convolution of two functions is defined as [Ref. 4 and 8]: 


mer) = r(t) * s(t) = I olson > , r(t-x)s(x)dx 


The cross-correlation is: 


Cc oO 


R(t) = Sf rx)s(x-t)dx = Ss r(x-t)s(x)dx 


OO =O 


This is a similar integral to the convolution with the difference 
that now none of the functions is time reversed. 
If one of the functions is taken, say r(t), and time reversed before 


calculating its convolution 


Be) = r(-t) * s(ty= sf r(xer}s(x)dx 
then 
g(t) = R(t) 


and it is concluded that the evaluation of both integrals reduces to a 
Single problem, since one of the functions can always be reversed prior 
to taking the product and integrating. 

The solution of these integrals is by no means a simple task. Even 
in the case where both functions are limited in time so that the integral 
vanishes outside a certain interval, that is, when it is only necessary 


to evaluate 


Sie 





i r(x )s (x-t ) dx 
“) 
the computation is not straight-forward. 

In practical applications this integral is normally evaluated by 
graphical or numerical methods. In both cases the approximate solution 
is still time consuming in general. 

Another method to evaluate the convolution integral can be reached 
through the theorems of Fourier transformation. 


Fourier transform and its inverse are defined as: 


Fw) = FLf(t)] = f f(x)eJ°%ax 


iH 


f(t) = FLFw@)] =o S Flw)el*tay 
The inverse transform of the product of the transforms of each of 
the functions r(t) ans s(t) is as follows: 


12 @) 
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2 fe ees r(xje Xx] S(w)erT dy 


oT 


oO 


a r(x) dx [a , S(w ele t-X) 43 = f r(x)s(t-x)dx 


00 


I 


Since, as is easily seen, the Fourier transform of r(-t) is R*(w), 


it immediately follows that 


fee) 


FoUER*(w)S(w)] =F r(x-r)s(x)dx 
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What these expressions denote is that if the Fourier transform of 
both functions is known, the evaluation of the convolution or correlation 
integrals may be obtained by multiplying the transforms and then performing 
the inverse transformation. This process greatly simplifies the task if 
there is a method to easily compute Fourier transforms. 

Such methods are not known, in general. However if the functions are 
represented in the discrete domain, there are fast solutions that can be 
obtained with digital machines. The computation then gives very close 
estimates for correlation and convolution. 

In the discrete case — the functions are represented as r'(t) and 


s'(t) thus [Ref. 5]: 


foe) 


r‘(t) = 2 r(nT) and s'(t) = 


N==-co 


s(nT) 


Bm g 


where T is the spacing between samples, and if the values of r(t) and 
s(t) are zero outside a certain interval, say 0 to (N-1)T as always 
happens in practical situations, 

N-] N-] 


r'(t) = &  r(nT) SEE ae es ine) 
n=0 n=0 


The discrete Fourier transform is defined as 


N-] 
R(K2) = = xr(nT 
n=0 


ye JeTnk 


with 2 = ce and an inverse transformation 


N-] 


E R(ka)el@tnk 
k=0 


r(nT) = v 
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The discrete Fourier transform has several properties, many of them 
paralleling those of the continuous transform. One fundamental distinc- 
tion, however, results from the definition of 2: the transform or the 
inverse transform has only N distinct values for k(or n) between 0 and 
N-l with these values repeating themselves for k outside this range. 
This may be seen as the transformation giving the result of a periodic 
sequence with period N. 

A similar relation exists between the convolution and the product 


of the transform of two functions as occurred in the continuous case: 


N-] 


f(nT) = 7 _ R(K2)S( kx ed 
N-] N-1 N=] , 
-7 Fee ore ets eetnnene peut 
k=O 2=0 m=0 
N-1] Ne] N-] : 
= 5 d, ye xr(2T)s(mT) [2 ed (n-m-£ )k24 
£=0 m=0 k=0 


It is not difficult to show that the summation inside the brackets 
is zero unless m is equal to (n-£) modulo N when the result is zero. 
Hence, substituting (n-2) modulo N for m: 
N-1 
f(nT) = £& ~~ r(&T)s((n-2)T) 
r=0 
where (n-2) is to be considered as (n-2) modulo N. f(nT) has then an 
expression similar to that given for f(t) and it may be said that f(nT) 
is a convolution. However it has the special characteristic of repeating 


itself with period N and so is named a circular or periodic convolution. 
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As a conclusion, an estimate of the convolution of the functions 
mt) and s(t) is f'(t) = a f(nT), where r and s are assumed periodic 
with period NT. - 

To evaluate that expression it is possible to take the discrete 
Fourier transform of each function, multiply them together and then 
compute the inverse transform. 

Fortunately, a very convenient and fast algorithm known as Fast- 
Fourier-Transform has been developed and is specially adapted for use 
with digital machines [Ref. 3]. 

Thus, it is now possible to evaluate convolutions and correlations 
for practical, real-time applications. 

But not all functions dealt with are periodic. On the contrary, 
very often they are limited to a relatively short interval of time. The 
method can still be applied if the functions are prepared before per- 
forming the convolution. Figure 2.4 illustrates how this can be done. 

The convolution of two rectangular pulses is shown in Fig. 2.4a and in 
2.4b the same thing is presented in the case of periodic repetition of 

the pulses. The periods were purposely chosen to space the pulse widely 
so that, when one series is slid across the other, there is no possibility 
of any single pulse of the first function coinciding with more than one 
pulse of the second function at the same time. One period of the convo- 
lution has then the same form as the convolution of the aperiodic functions. 
Therefore, in order to evaluate the aperiodic convolution it is necessary 
only to extend both functions with zeros, constructing periodic functions 
with a period size of MN, where M and N are the actual sample sizes of 
the original functions, and taking only one period from the resulting 


periodic convolution. 
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CONVOLUTION 


a) Aperiodic 


b) Periodic 


FIGURE 2.4 
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Another application of this method refers to the case when one of 
the functions is very short in comparison with the other, as appears in 
a detection problem when a short replica is correlated against a very 
long incoming signal. This situation is presented in Fig. 2.5. 

The long sequence is partitioned in sections of size M, each section 
extended with N points from the previous section, and the samll one 
extended with M-N zeros. If the correlation of the N samples (plus M 
zeros) against successive N+M sized sections is performed, and the results 
added, it is evident that the first N points of each partial result will 
be originated in the same way as the last N points of the previous sec- 
tion. This distortion is avoided if the N first points of every partial 
result are discarded. 

In the derivations mentioned above, no restriction was imposed upon 
the functions that can in general be complex. Therefore in practice a 
group of four convolutions or correlations is performed at the same 


time, 1.@€., 
r(t) = ri(t) + dr (t); s(t) = s.(t) + js, (t) 
x _ * ae * * ; * 
moeits(t) Gs. 7 Na Se Tdi Se Tees 
In most applications, however, real functions are dealt with and it 
is then possible to perform more than one operation simultaneously. For 


instance, the partial result can be evaluated for two sections in one 


Operation, or: 


r(t) = sae j0 s,(t) = 5s. + j0 s,(t) See j0 


Ir 
= : * = * iy x 
2c 1p 7 J5or r(t) s(t) "ry lr ‘i JT > 2p 


After the operation is concluded the real part is taken as one section 


of the result and the imaginary part as the other. 
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The subroutine FASTCOR uses for the evaluation of the Fourier trans- 
form a very fast program, FOURZ2, based on the Cooley-Tuckey method, which 
is available in the library of the Computer Laboratory. This program 
requires that the functions to be transformed have a number of points 
equal to a power of 2, and is the reason for this restriction on the 
sample size. 

FASTCOR has three main options either for convolution or correlation. 
They are assigned by the argument IT, i.e., periodic, aperiodic, or 
recursive. Figure 2.6 shows a flowchart for FASTCOR. 

Initially the size of the samples is checked against the option 
selected. In any case, the prepared sample (a ascent) so: both samples 
has to have the same number of points. 

If a periodic operation is chosen and the second function is longer 
than the first one, with a maximum of 512, it is reduced to the shorter 
length. If an aperiodic operation was selected and the sizes of both 
Samples are equal, in which case the maximum length is 256 (M), the two 
functions are padded with zeros to 2M. In case the first sample is 
Smaller, it is limited to one-half of its size and the second one has 
the number of points reduced to a value equal to its original value 
minus one-half of the first function's original size. Both functions 
are then padded with zeros up to the initial size of the second function; 
this assures no overlapping in the result. Finally, with the recursive 
option the first function is always assumed to be smaller than the second 
which has a maximum acceptable length of 512; the first sample is padded 
with zeros to the length of the second. 

The next step is to flip the first function, send it to be displayed 


on the screen, since it will soon be erased from memory, and store it as 
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complex with imaginary part zero. If a correlation is to be performed, 
the function is displayed without reversal but reversed and stored as 
complex; for convolution the opposite occurs. 

If the option is not recursive, the second function is also stored 
as complex with zero imaginary part, and both are then Fourier trans- 
formed, multiplied together, and the result inverse transformed. During 
these operations some transfers of storage occur to minimize the space 
used. Finally, the second function and the result are shown on the 
screen unless no display is required. 

When the option is recursive, an auto-correlation is initially per- 
formed on the first function to establish the normalization coefficient 
and the actual threshold value. This normalization factor is computed, 
Since the maximum amplitude a correlation can attain occurs when the 
two functions are identical and therefore corresponds to the maximum 
value of the auto-correlation. But since the amplitudes of the two 
functions have been previously and independently normalized, the new 
normalization factor is corrected by the ratio Q mentioned earlier. It 
was chosen to normalize the threshold value to the maximum amplitude of 
the auto-correlation, which means that a relative threshold of one would 
correspond to a perfect matching of the functions. This actual value of 
threshold is then further normalized, in the same way as the result, for 
display purposes. 

Having completed this task, the second function is then stored as 
complex. In the strictly recursive case a double-sized (2M) sample is 
stored as real, and the next section, with overlap, as imaginary. For 
the determination of probability of detection, however, a sample section 


padded with zeros is stored as real, and the next section as imaginary 
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so that effectively a double aperiodic correlation is taken at one 
time. After correlation, the two sections are moved in core, compared 
against the threshold, and displayed if so desired. 

If a correlation is being evaluated, the threshold level is also 
displayed, and an indication appears as to whether this amplitude is 
reached or not. Information also comes to the screen about the probabil- 
ity of detection that has meaning only if this option is in execution. 
The last measured input signal-to-noise ratio is also shown. 

Every time the threshold is met in a distance measuring operation, 
the actual distance, represented by the sample number, is displayed. 

The signal, in this case, is added to noise only at one predetermined 
distance, 1.e€., only in one group of samples. 

After this phase is completed the first 2M samples of the second 
function (signal) are discarded and the last M are translated to consti- 
tute the new first set. SEGSAMP is then called with INDEX equal to 3 
to fill out the function with 2M new samples, and the computation starts 
again. 

Exit from the recursion loop results when either the loop has been 
executed a preset number of times or when the operator presses a switch 
at the graphics console. 

In all cases the control returns automatically to the sequence 
indicated by LABEL except when the probability of detection (or the false 
alarm rate) was requested in the ANALYSIS mode. In that case the proba- 
bility of detection (or the false alarm rate), threshold setting, number 
of sections processed and input signal-to-noise ratio are all printed 
on the line printer, and the operator must command continuation by 


pressing the RETURN key of the terminal keyboard. 
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F. SUBROUTINES ANALYSIS AND SEGFOUR 

These subroutines belong to the ANALYSIS mode. 

ANALYSIS is intended to clarify further the request of the operator 
and route the computer to the appropriate sequence. The argument IENT 
is a pointer to indicate whether the subroutine is entered to initiate 
a sequence or to end it. This is a requirement of the overlay structure 
Since the space it occupies in core has to be released whenever FASTCOR 
is called. 

The initial selection is indicated by IJ, kept in common storage. 

IJ equal to one indicates that sampling is desired. Before calling 
SEGSAMP a reminder is displayed asking the operator if the constants 
are adjusted to sample the correct function, for the correct processor 
and from the correct point in the analog (INDEX, ISAMP, LINE). If the 
answer is negative the constants are displayed to allow corrections 
after which sample has to be requested again. An affirmative answer 
initiates the sampling operation and then the new sample is saved in 
background storage. 

IJ equal to two indicates Fast-Fourier-Transform and the operator 
is requested first to indicate which function: replica or M sized 
function, signal or 2M sized function, correlation or the result of the 
last operation performed by FASTCOR. Then SEGFOUR is called to prepare 
the function selected, perform the transformation through FOUR2, calcu- 
late the absolute value and the phase angle of the transform, send them 
and the original function to the display and read back from’secondary 
Storage the functions destroyed during the operation. 

When correlation is selected, IJ=3, it is only necessary to indicate 


if the mode desired is periodic, aperiodic or recursive. FASTCOR is now 
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called through FASTBRG, a bridge subroutine required to save addresses 

and parameters while the segment swapping takes place. Back from FASTCOR, 
FASTBRG re-enters ANALYSIS to restore the functions destroyed and put the 
result of the correlation in backstorage. 

For IJ equal to four (convolution), five (power spectrum), or seven 
(auto-correlation), the operator must select the function and mode. With 
convolution, the “other function" is always assumed to be the replica; 
with the other two operations, recursive is an illegal choice and if 
made results in an error message. The functions are then moved in 
storage so as to send to FASTCOR the appropriate parameters. Returning 
to ANALYSIS, action similar to that for correlation is taken, except when 
power spectrum is the operation. 

The evaluation of power spectrum is done through the application of 
the Wiener-Kinchin relations [Ref. 8]. 

It was earlier recalled that the convolution is the inverse Fourier 
transform of the product of the transform of the functions to be convolved 
and that the correlation can be treated as the convolution taken with 
one of the functions time-reversed prior to the operation. Also the 
transform of r(-t) is equal to the complex conjugate of R(w). Now, 
the auto-correlation of a signal for which the energy content is not 


finite is defined as 


P(t) = at u p r(t)r(t-t)dt 
ie 
2 
Hence, 
P(r) = 22 Atr(t)*r(-t)] 
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where G(f) is the power spectral density. If the energy content is 
finite, the energy spectrum is considered instead of the power spectral 
density. Therefore, since the power (or energy) spectrum and the auto- 
correlation form a transform pair, the spectrum may be evaluated by 
computing the auto-correlation and taking its Fourier transform. 

ANALYSIS then determines the auto-correlation to be performed and 
when re-entered, calls SEGFOUR to transform the result. (None of these 
results are saved.) 

The last control of the subroutine corresponds to IJ equal to six 
or the display selection. After the answer for the function is selected 
and decoded the subroutine DISPLAY is called to show that function. 

ANALYSIS has to keep track of the number of samples related to each 
selection and controls the allowed limit for all situations reducing the 
Size of the sample if necessary. The common location II is used for 
this purpose during the overlays substitution. The size of the last 
result (CORR) stored is saved in IMEM. Its value is set when the sub- 
routine is re-entered and corresponds to the number of points displayed 
as indicated by LBG(3). 


The flowchart for this subroutine is shown in Fig. 2.7. 
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G. SUBROUTINES PROCESS AND PISPLAY 

PROCESS is constructed to chain various sequences in order to execute 
the operations necessary to compute input signal-to-noise ratio, evaluate 
the processor input-output performance, and determine the ROC curves. It 
also permits simulating a situation of target detection and distance 
measurement. This automatic procedure constitutes a part of the SIMULATION 
mode except for the signal-to-noise calculation that may be called 
separately in the ANALYSIS mode. 

Except for distance measurement, the program establishes loops where 
some parameters are varied for every recursion until a limit is reached. 
Programwise, however, these loops are interrupted several times by trans- 
fer to another overlay. More common storage is then required to save the 
Status during the segment swapping, and to allow resumption of the loop 
an the proper position. 

A simplified flowchart for these sequences is presented in Fig. 2.8. 

When signal-to-noise ratio is called in the ANALYSIS mode, it is 
initially necessary to scan the signal and noise injection potentiometers, 
Save their settings, and restore them after the computation is completed. 
Entered the first time for performance or ROC determination, it determines 
the number of steps that will be used to calculate the SNR. The operation 
really starts by saving the processor status and establishing a linear 
processor, since the input signal-to-noise ratio (after the transducer) 
is desired; the processor is reset after the computation. Two loop paths 
are then executed, one to determine the signal power at that level of 
injection (no noise present), and the other for noise power (no signal 
present). The ratio, in decibels, is taken and stored as SNR, and 


displayed if the operation CONST is selected next. 
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The Wiener-Kinchin relation again: 


Piz) = f£ G(f) etc™ ‘ae 


When t iS zero 


P(O) = i G(f) df 


oO 


meaning that the value of the autocorrelation when t = O is the total 
average power. 


The absolute value of the autocorrelation is 


P(t)} = [ey @ctyede™ Tae] < 5 a(teJe™ ae] = sy G(F)df = P(O). 


Pao 0) ee] 


So, the power content can be found by taking the maximum value of the 
autocorrelation since P(0) is an absolute maximum. This is exactly what 
the program does. 

The input-output performance or processor gain can be defined by 
the equation 


SNR, = f(p,W)SNR, 


where SNR, is peak output signal power to output noise power ratio 
SNR. is the input signal-to-noise power ratio, and f(p,W) is a function 
of the processor and its associated information bandwidth. 

The performance operation begins sampling and storing the replica 
in a position never to be destroyed during the process in this particular 
case. Then it uses the previously established values for signal and 


noise injection that correspond to the highest input signal-to-noise 
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ratio to be considered, and sequentially computes input signal-to-noise 
ratio, samples noise alone, cross-correlates replica and noise, auto- 
correlates this output and takes the maximum (output noise power). It 
then samples signal plus noise, cross-correlates replica and signal plus 
noise, finds the maximum (signal peak output). Next it reduces the 
Signal injection and correspondingly increases the noise by a pre- 
determined step and recycles as long as the signal-to-noise ratio is 
greater than zero. During the operation the individual values of signal- 
to-noise ratio, output noise power, and signal peak-output are placed in 
common storage. When the loop ends, the computed values are rearranged 
in core and SNR, is computed. Finally, the subroutine PISPLAY is called 
to display the resulting curve with the abscissa as SNR,(db) and the 
ordinate as SNR, (db). (The input signal-to-noise ratio is measured after 
the transducer filter. ) 

One detail must be examined in the output noise power computation. 
Since the cross-correlation of replica and noise has 2M points, it may 
exceed the maximum limit allowable for autocorrelation if M is 256. The 
solution adopted was to calculate the autocorrelation in two steps, 
one-half at a time and then add the two maxima; this is legitimate since 
tne correlation is a linear operation. 

For ROC determination, the procedure, up to the input signal-to-noise 
ratio evaluation and the recursion for the several levels of signal, is 
identical to the sequence just described, with the exception that the 
replica is placed in secondary storage. Inside the loop however two other 
loops exist: the outer one for the threshold value that steps down in 
0.1 steps from the initial setting until zero is reached; the innermost 


loop has two paths, one for the calculation of the probability of 
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detection (signal plus noise) and the other for false alarm rate (noise 
only). PISPLAY is also used to present the results on the screen. The 
abscissa is presented as false alarm probability and the ordinate is 
false alarm rate (per sample). A listing of the signal-to-noise ratio 
for the curves presented and normalized threshold are also displayed. 

According to this description, instead of the receiver operating 
characteristic the result is a modified ROC where false alarm rate was 
substituted for false alarm probability. This modification has more 
practical significance and is also easy to evaluate in a short period 
of time while the probability of false alarm is not. These two param- 
eters are related by the noise output bandwidth. | 

At this point a warning should be made. The above estimations serve 
as an indication of the probable behavior of a processor and that is the 
purpose of the system. More pDrecise ensemble average estimations would 
require lengthy repetitions, and it was not considered worthwhile to 
incorporate them. However, as will be pointed out later, the basic idea 
can be used for a much less flexible fast analytical system or for 
Special purpose designs. 

Detection and distance measurement form the last option of PROCESS 
which sets parameters and lines to initiate the procedure, and samples 
the replica and the initial sections of signal. A line is set to dis- 
connect the pattern generator. ICOUNT is set to allow a maximum distance 
Scan corresponding to 20M samples and ISTF to connect the signal at a 


distance of 12M sampling intervals. 


H. SUBROUTINE CIRCPAT 


The SIMULATION mode is completed with this subroutine which is able 


to respond to the typed commands and set circuits and parameters to 


ie 





Simulate different transmission signals, such as CW burst, FM slide, 
FSK and propeller signature. 
These simulations are better discussed in the section dealing with 


applications. 


I, AUXILIARY SUBROUTINES 

A number of auxiliary subroutines were developed to enhance the 
utility and convenience of this system for the investigation. 

GRAPHAL has already been referred to and permits the control of the 
analog computer from the graphics terminal. It is entered whenever the 
command for one of its functions is typed; if this function is to scan 
the analog, ANOSCAN is called. These subroutines may be modified for 
general uSe. 

PPLOT copies the functions displayed on the screen and has three 
options to output them. The first is through the Tibrary routine VPLOT, 
with identifier changed to WPLOT since it was necessary to have it loaded 
only in the overlay in order to reduce core storage requirements. It 
plots the curves on the line printer in a one-page sized plot where the 
X-axis is transverse to the paper. NONGPLOT, the new identifier given 
to the library subroutine LONGPLOT for the same reason as above, also 
prints on the line printer but with the x-axis longitudinal to the paper 
and a plot occupies a variable number of pages depending on the number 
of points plotted and on a parameter set by the operator. Through another 
sequence, PPLOT uses the eight channel recorder as the output medium. 

When selected, it first asks that the recorder be made ready, and when the 


operator presses the RETURN key it waits 20 seconds before starting the 
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plot. The information about the potentiometer PO1i displayed there 
indicates the normalized threshold level that, if desired, may be intro- 
duced in that potentiometer. For this case, an analog probability counter 
will deliver pulses to the recorder whenever the correlation curve crosses 
the threshold value. The value for adjustment of POI1, however, is 
meaningful only if the curve displayed corresponds to a recursive type 

of correlation. 

Subroutine DEMANAL routes the command to other subroutines which 
perform correlations through a direct approach, thet is, by simple inter- 
pretation of the integral definition. Those subroutines are: ANACOR 
which performs a periodic autocorrelation using an analog multiplier and 
integrator; DIGCOR which evaluates auto or cross-correlation entirely 
digitally; and RECOR which computes recursive cross-correlation. 

Subroutine NORMALIZE is intended to normalize a function to some 
preset value. The storage address, the maximum value desired and the 
number of points are the arguments passed to it. 

When the program is first loaded or after GRAPHAL is called, INITIAL 
initializes the text and graph directory and the graphics terminal. It 
is entered to erase the text blocks after an operation is concluded or 
to erase the graph blocks after a new selection is made (except if the 


selection is CONSTANTS or any of the plot options). 


J. SUBROUTINES DECO AND CONST 

DECO works as the control section of the system, allowing the user 
to type in his selection which is then decoded and interpreted to channel 
the actions in the desired direction. It forms a table of control words, 


composed of the first four letters of the selections available, and when 
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a command is decoded checks it against that table and sets two pointers, 
II and IJ, and returns to MAIN. II is then used to select the segment 
while IJ indicates the correct choice. If no entry exists in the table 
matching the selection, an error message is displayed and a new order 
may be entered. DECO is called automatically after the completion of 
an operation or when an illegal procedure is requested. in another sub- 
routine. 

The control function is completed with CONST which displays various 
parameters and pointers and permits the operator to alter them inside 
allowable limits. The names under which they appear are not necessarily 
those used in the program, but others more easily eenenbened: A full 
discussion of this will be given with the graphics terminal usage 


explanation. 


is 





IIT. ANALOG CIRCUITS 


ine analog computer was selected to simulate the signals because 
they are generated in the same form as they occur in reality. Thus the 
analysis of existing signals can be made by the simple substitution of 
the actual signals. In this case the analog computer itself acts as 
the interface to the system. 

Another important feature considered was the relatively large 
parallel logic capability of the COMCOR Ci-5000 analog computer, which 
allowed construction of control circuits to enlarge the range of utili- 
zation and analog processing. 

The analog computer circuits included: one sine generator used as 
PATTERN GENERATOR, one HARMONIC GENERATOR fed by three independent sine 
generators, one WHITE-NOISE generator, a SIGNAL COMBINER, a BANDPASS 
FILTER, an INTERRUPT GENERATOR and some auxiliary circuits. 

As a general rule the circuits and their parameters are remotely 
controlled from the graphics terminal with exception of those few which 
need patch relocation. All lines that can be set are tied to another 


test line in order to permit status check from the display. 


A. PATTERN GENERATOR 


1. Theory 
. 2 
The differential equation y = oy = sas} represents a sine wave 


dt 
and its solution can easily be implemented with integrators and adders. 


The mechanization of this equation is presented in the following expres- 
sions where the bar means that the variable is a scaled (computer ) 
variable, theg's are amplitude scaling coefficients, 8 is a time scaling 


coefficient, and t. is comoputer scaled time interval: 
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2. Implementation 


The basic implementation of the equations above is shown in 
mig. 3.1. . 

Since it was desired to make this circuit simulate a CW burst, 
as well as a linear FM slide and an FSK codes, a multiplier was utilized 
for P and Po because it allows changing frequency during operation 
without resetting the amplifiers. Additional circuitry was required to 
Generate the multiplier reference voltage, to select the type of opera- 
tion, and to switch the generator on or off. The final diagram is 
presented in Fig. 3.2. AQ1l7 is the integrator for y and A021 for -y. 
A203 and A204 provide polarity reversal as required. The multiplier 
M401 implements P, and M402 introduces Poe Relay R26 connects the 
reference input of the multiplier either to the ramp generator A025 or 
to the FSK relay R10. P030 and P031 adjust the lower and upper shift 
reference using A205 and A206 as loads. R006 is the relay that switches 
the generator output on or off. 

The ramp generator computes the value of R = I + St where R is 
the reference voltage, I the initial setting and S the slope. P026 
establishes I, P025 fixes S and P027 fixes the threshold value for CO/ 
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which controls the reset time for the ramp. P032 and C04 are responsible 
for the shift instant. P02] prevents the necessity to set P025 to very 
small settings and A027 and A030 serve as load to A025 when R26 is open. 

To make it possible to run this generator at very high actual 
frequencies (150 Hz), positive feedback is introduced by P017. The 
setting of this potentiometer is not critical since the limiter LOO 
prevents the amplitude from growing. 

3. Logic 

Fhe logic circuitry required to control the pattern generator, 
presented in Fig. 3.3, is the result of the logic equations developed 
in the Appendix. 

Lines TOQQO7 and TO10 are set to control the integrator mode 
through several drivers. 

The NAND gates, comparator C07, delay-flop DFO2 and line TOl1 
direct the ramp generator operation. It then follows the general mode 
except if CW is required when it remains always in reset. During compu- 
tation the ramp generator is made to reset every time a sample period 
ends. 

C04 actuates the FSK switch R10. Lines T1016 and TO17 bring the 
voltage to R26 and RO6, respectively. 

4. Scaling 

The presentation following has the purpose of establishing simple 
reference tables for quick problem scaling. The final objective is to 
find the value of R, the reference voltage, and the amplifier gain suit- 
able for each case. 

The amplifier gain (G) is a function of the input conductance 


(I-mmhos) and the feedback capacitor (C-nfd), coefficients p, and p, 
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mentioned earlier are the product RG divided by the constant 100, 
characteristic of the multiplier. For convenience, R was multiplied 
by 100 in the final result to become the integer value that is sent to 


set the potentiometers. 


where a and @ are amplitude scaling factors for y and y, 6 the time 
scale factor, ts the computer scaled time interval and ire and We the 
scaled frequency and angular frequency respectively. 

Since y = A sin wt, y = -Aw cos wt assuming A equal 100 for a 


maximum voltage scaling a = 1 and a =w. Thus 


_ 100 x 100 C w _ 20000nf C 
BI I 


3 


R has to be kept below 10000 (100v) and it was chosen to keep 


R 


C as .| for real time and .01 for time scale for an actual computer 
frequency, Ts around 15 Hz and 150 Hz respectively. The other values 
follow as a consequence and are tabulated in Table I in the Appendix. 
9. Operation 

For CW burst code, A027 is kept in RESET through 1011, R26 is 
held in the normally closed position (NC); P025 is set to zero and P026 
to the reference voltage R. 

In FM slide A027 is made to compute once every sample and then 
reset. The relay R26 is maintained in the normal ly closed position(NC). 
The potentiometers are set in the following way: 


P02] to the inverse of the gain of A025; 


P025 to the ratio of the upper reference voltage(Ry) minus the lower 
reference(R, ) to the actual sampling time, j.e., 
(Ry-R, ) x sampling period x number of points in the sample. 
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PO026 to the lower reference voltage; and 
PQ27 to the upper reference voltage. 


: For FSK the same adjustments as for FM slide are necessary with 
the following alterations: 


R26 - placed in normally open (NO) through the action of T1016; 
PO30 - to the lower reference voltage; 

P03] - to the upper reference voltage; 

PO32 - to the central value of the reference voltage R. 


In all cases above, if distance measurement is required the 
relay RO6 is maintained in normally open (NO) all the time except at the 


moment when the burst is to occur when it closes under control of TO17. 


B. HARMONIC GENERATOR 
me sine generators 
The three sinusoidal oscillators included in the harmonic generator 
are of the basic type presented in Fig. 3.1, with P and Po adjusted by 
potentiometers. The scaling formulas for the multiplier case also apply 
ener 
The harmonic generator was built to simulate a relatively low- 
pass line spectra similar to the signal generated by a propeller. One of 
the sine generators produces the signature while the others produce a 
"non-white" noise, as if they were other identical sources of slightly 
different frequency. 
Zee Mul tipliers 
In the diagram of Fig. 3.4 it is seen that the three generators 
were combined and fed to four multipliers arranged to give as output, 


after combination in A052: 


7 2 3 4 5 6 
E = Kye + Kee + K€ + Kae +: Kre + Kee 
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The parameters K are adjusted in order by P051 to P056 with the first 

three feeding the amplifier at I = 10 and the last at I = 1. This 

Peauires that POQ51 through P053 have a maximum of 10 while the sum of 

all three, times 10, plus the remaining three be 10000 or less. Any 

one of the parameters may have its polarity reversed by the switches 

ROO, ROZ, RO4, R20, R22, R24 activated by the lines TOO] through TO06. 
Table II in the Appendix gives some of the combinations of coef- 

ficients and frequencies up to sixth harmonic when one or two frequencies 


constitute the input signal. 


C, WHITE NOISE GENERATOR 


1. Theory 


To build a source of random noise a pseudo-random sequence 
generator was used. 

As is well-known, a true random binary sequence may be classified 
by three properties [Ref.2]: 


a) The number of ones approximately equal to number of zeros; 


b) About half of the no-change states have the length of one, 
about one-fourth have length two, one-eighth have length three, 
and in general 1/2" have the length n; 


c) The expected value of the autocorrelation has a peak at the origin 
and decreases rapidly away from there. 
A linear maximal sequence, like those generated by specially 


designed feedback shift registers, has properties very close to these: 


a) The number of ones is exactly one more than the number of zeros; 
b) The number of the no-change states of length n is 1/2", 
c) The finite autocorrelation, the sequence considered periodic, 
has the form of 
L 


P(k) = = aa , » where L is the sequence length and a. the 
value of the nha point of the sequence. This expression has only 
two values, say L when k = 0, +L, + 2L, ... and -1 otherwise. 
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Thus, within appropriate limits, a pseudo-random sequence may 
be used to simulate a true random one. 

To generate binary analog pseudo-random noise from a random 
sequence, two techniques can be used. In one of them all bits are added 
modulo 2 for each shift. For the other only one bit is taken in each 
shift. This logic signal is then amplified and clipped to originate a 
pseudo-random square wave. The second technique was adopted. 

The analog noise generated has first and second order statistics 
essentially identical to those of true random noise if the time of 
observation is short compared with the sequence period. Also, it is 
found in the literature that, after lowpass filtering, high confidence 
level tests indicate essentially a first-order Gaussian behavior [Ref. 7]. 

Again, for a long sequence and a relatively short observation 
interval, it is not possible to resolve the spectral lines that neces- 
sarily exist in the power spectrum of a periodic waveform. 

2. Implementation 

As shown in Fig. 3.5, a 25-bit feedback shift register (FSR) was 
constructed with bits 3 and 25 being added modulo 2 by the quarter-adder 
formed by A400, A200, A100, and then fed back to the input. Tnis arrange- 


cD 1 = 33, 554, 431. 


ment generates a maximal sequence with length L = 2 
The output of the FSR and its inverse were used to actuate 
digital-analog switches, DAO] and DAQ2, which apply -100v or +100v, 
respectively, to the amplifiers AQ1O and AQ16, whose outputs are further 
combined in AQ14. 
In order to avoid the all-zero state, bit one is set at reset 
time and all the remaining reset. 


One requirement of a system built for signal analysis is the 


possibility to return to any observed situation under operator control. 
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For this reason, the shift pulse and the reset lines were made logically 
controlled by a special circuit, to be described later, so that the FSR 


operation is synchronized to the generator modes. 


D. SIGNAL COMBINER 

All three types of signals just described were brought to an 
amplifier, AQ24, through the potentiometers, PQ22, P0023, P0024, as shown 
in Fig. 3.6. The amount of injection of each one can be controlled and 
any combination achieved. 

Potentiometer PQI5 was included to feed an initial condition to AQI5, 
if desired. With this feature and no inputs to AOI5 and A024, it is 
possible to simulate an impulse to the filters that follow and then 
analyze their characteristics. 

Line T1500 was connected at the output of the combiner to make it 


possible to sample before the input filter if desired. 


E. BANDPASS FILTER 

Actually, two filters were included after the combiner to simulate 
the transducer and the input circuitry. The original intention was to 
have a narrow-bandpass filter for the active detection simulation and 
the lowpass filter for the propeller signature problem. 

1. Narrow-band Filter 

A second-order filter, shown in Fig. 3.7 was installed for this 
purpose. Its frequency response is given by 
jwA 
: (ww) + jwB 


Ww) 
where A is the gain and B the bandwidth related to the Q by B= 7: 
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The differential equation describing this circuit is: 
oe - 2 i é 
y + ABy + Aw, y = -Ax 
which, when integrated and put in terms of computer variables, becomes: 


yas 


ToD [to 


Pt 2 dt. - f= x dt 


WW 
Z= ~fe (-y) dt, 

PQ36 then serves to adjust the gain, P037 the bandwidth, and 
PQ40 and P04] the central frequency. The values of Table I may also 
be used for these potentiometers if the frequency is substituted by 
BA/2n, Af and Fo respectively. To avoid small settings for P036 this 
potentiometer was normally connected to an amplifier having a gain ten 
times smaller than that used for the fy input. 

The filter characteristic may be estimated by application of an 
impulse, as just mentioned. An example of this was taken and is shown 
imerig. 1.9. 

Line T501, connected at the filter output, makes it possible 
to take samples at this point. 

2. lowpass Filter 

A simple RC filter was used for lowpass. 

To assure the first-order Gaussian characteristic, the time 
constant RC has to be much greater than one shift period, At, but still 
less than 25At, the sequence period. Of course, the inverse of RC is 
the cutoff filter frequency. 

Table III in the Appendix was constructed to facilitate the 


adjustment of R, C, and At or the clock pulse frequency. To extend the 


2) 








range of possible clock frequencies and make it easier to meet problem 
conditions, the built-in computer clock was passed through a frequency 
divider with outputs equal to one, one-half, one-fourth and one-eighth 
times the input. 

If the narrow-band filter is in use, the shift pulse is set to 
make At much smaller, and 25At greater, than the inverse of the upper 


frequency. 


F. INTERRUPT GENERATOR 

This was one part of the system treated with very close attention. 
Since the digital computer program may be destroyed by interrupts arriving 
at the wrong times, they had to be carefully inhibited except when they 
should be expected. 

Figure 3.8 shows the interrupt generator. As a general rule, the 
pulse that causes the internere is placed at the input of a NAND gate 
that has also an input connected to a line that is kept off by the 
program except during the sequence when it is legal for the occurrence 
of that particular interrupt. In this situation the output of the gate, 
if a pulse arrives, actuates a delay-flop that effectively generates the 
interrupt. Delay-flops were used to assure a short, well-shaped pulse. 

Another interrupt generator is the circuit formed by the lines TOI3 
(set), T113 (test) and the switch to ground. In the subroutine FASTCOR 
the digital computer tests line T113 once in every recursion. If the 
operator closes the switch at the display console, the line is set low 
and the program, after reading this value, branches out of the recursion 


loop, which corresponds to a software generated interrupt. 
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G. AUXILIARY CIRCUITS 

In addition to the functional circuits described above, logic 
circuits were developed to insure the proper operation of the total 
system. 

Among these, the synchronizer (Fig. 3.9) is responsible for keeping 
the noise generator in the same mode as the other generators. This 
circuit also prevents the occurrence of sample interrupts when the 
generators are not in a COMPUTE made. In the Appendix are included the 
logic equations represented by the synchronizer. Other circuits and 
connections are shown in Figs. 3.10 and 3.11. 

The real-time/time-scale controls were placed under digital command 
through the lines T030 and T031, respectively. 

A four-bit counter was patched to work as a divider for the built-in 
clock to facilitate the rate selection for the feedback-shift register. 

As a sample clock, a frequency synthesizer with square wave output 
was used. An analog comparator fed by this square wave generates a 
pulse starting at the zero crossing to trigger the sample interrupt, or 
to function as a timer for the recorder copy function, RPLOT. 

The eight-channel recorder was connected to the outputs of the 
Signal combiner, two of the sine generators and the pattern generator. 
Channels 6, 7, and 8 receive the copied curve through the DAC lines. 
Channel 5 displays a pulse whenever the third curve being plotted reaches 
the threshold setting; this pulse is generated by the probability counter. 

The last circuit built was a multiplier-integrator for secondary 


computation of correlations by the subroutine ANACOR. 


94 















A0GOW 


Ol4G 
cee oe 


345 9SNa 





OlIN 


INIT LSY 3Q 19 


(AIdWVS) OPN OL 
lOZN 


SNIT A419 


Z.0IN 


Ore ails 


IOYLNOD 
JGOW SOIVNY 





OIZN 


SOON 


YSAZINOYHONAS SAQOW 


Q34u 


vOodWd 
G00 W4 


Y3AQIAIA 
AION WS 


oI 





L2vil- €3AYND 
9Zbl - Z23ANND 
GZ2vl - 1 3ANND 

bvOV-LNNOD sOXd 
LIOV — N39 LLVd 
JOOV- N39 ANIS - ¢ 


F0OO0V- NAD SNIS -¢ 
pcOV — YSNISWOO - | 


j 
¢ MOK @O 


Y3CYOOSY 


OIPN P¥e 
SOCN OL 








MIO19 
TWWNYSsLxX4 


909 
AMIO1D AIMNVS 


Ol'¢ 3YNSIS 


Sdd 00! 


| Sddm | 
Eclaieeis|ieisajOclsaj \ 


$dom Ol 
SdOxA00! 


M90 19 
LON OL NI LA 


YSGIAIQ ADNSNDSAYS 


a} D 
wl ep. CD 


IOYLNOO 
31V9S SWIL/SWIL Wad 


ATO 





BAL MAL Aly 








We B3YyNSISs 


S 


GOO 


io 


¢ 


YSINANOD ALINIEGVEOud 


YOLVYOSLNI 


YSIIdILINW 


ay, 





TV. GRAPHICS TERMINAL 


The graphics computer provides the system one of its most important 
features, that of human-machine interaction. 

The real-time operation enhanced by the natural form of graphical 
presentation made it easy for the operator to assimilate information, 
make decisions, and select the next sequence using the system control 
capability provided by the graphics console, | 

Since the AGT-10 is a general purpose computer it may be used to 
alter either the data or the form of the presentation to optimize the 
man-machine interface. 

The system was controlled from the graphics terminal by typing simple 
messages on the screen using the console keyboard. A number of control 
possibilities are available. The operator could select to put the 
analog computer in any mode, or test its line values, select signals, 
sample, take the Fourier transform, perform auto- or cross-correlation 
and convolution, or determine power spectrum. He could also either 
initiate sequences to evaluate processing gain, ROC, distance, or signal- 
to-noise ratio or he may manually control every step of these operations. 
Finally, he had a choice to copy what was presented on the screen either 
on the line printer or on the paper recorder. 

In the Appendix is given a more detailed description of the operating 


procedures, 
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V. SONAR SYSTEM APPLICATIONS 


As a verification of the system capabilities several problems have 
been run to compare performances of the three processors, check the 
effects of integration time and bandwidth, and to understand the signi- 


ficance of information and joint information bandpass. 


A. PRELIMINARY CONSIDERATIONS | 


To reproduce a signal from its sampled representation, in an ideal 
case, it is necessary to take samples at least at the Nyquist rate of 
two times the highest frequency contained in the signal. In the case 
of a narrow-band-limited signal, when the bandwidth is much smaller than 
the highest frequency, the rate becomes twice the bandwidth. The 
recovery of the signal, however, requires filtering with steep filters 
not found in practice. 

For the system described, the only filtering action to which the 
Samples are submitted is represented by the manner used to generate the 
vectors for display. This action can be compared with that of a first- 
order interpolator. For this reason, whenever the shape of the functions 
was desired, a sampling rate was used of several times the highest fre- 
quency, even in the case of bandpass limited signals since no bandpass 
filter was implemented. 

For the time-limited and bandpass-limited functions the number of 


independent samples (n) is twice the time-bandwidth (BT) product as n=2BT. 


'Most of the expressions mentioned in this section are from, or follow 
directly from, those presented in the theoretical treatment of this subject 


by A. A. Gerlach, The Theory and Applications of Statistical Wave-Period 
MaOceSsmmonscook Electric, 1967. 
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At least this value is required to exactly represent the function and 
half of that if only the energy content is the object of the calcula- 
tions. When computing statistical measures of sampled functions these 
numbers must be taken into consideration; in particular, due to the 
relationship between energy and autocorrelation, it is expected that the 
number of independent samples be m= BT for this case. 

The practical limitation of the bandwidth B requires that when the 
Signal is processed no noticeable differences result if the bandwidth 
1s extended beyond that limit. This of course depends on the capacity 
of the processor to "use" or "recognize" the information conveyed. In 
the case of a matched filter or linear correlator this optimum bandwidth 


approaches the information bandpass W. [Ref. 4] where: 


oO 


rr o@(f)aty° 


=O 


CO 


fr 6(#)°aF 


=O 


In the general case of detection, not only the desired signal is 
present, but also the undesired noise, and the purpose is then to identify 
Signal in the background of noise. The effect of noise is to destroy part 
of the identifiable information or reduce the information bandwidth. This 
action is taken into account when the significant bandwidth is taken as 


the joint information bandpass or 


f G(f)df s N(f)dF 


Sif co 
S G(f) N(F) df 


OO 


From Shannon's expression for the information content, 1.e., 


m = BT log. (1+SNR) = [B log, (1+SNR) ub 
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this interpretation can clearly be verified. Therefore, the information 
content of a Signal in the presence of noise as seen by a cross- 


correlator type detector may be taken as: 


m = WT “where W is the joint information bandpass for that 
particular processor, signal and noise, and the total 
number of independent samples n= 2WT. 


It has been determined that for both the linear and the amplitude- 
limited cross-correlator, W= W, , 3 and for the clipper-1 imited 
processor W = 1.5 Neon In any case the computation of Ms ae is made 
with the signals as they are actually presented to the correlator, that 


is, undistorted, limited and clipped respectively. 


B. PERFORMANCE COMPARISON 
To compare the performance of the three processors, an FM slide 
Signal was chosen with the following characteristics: 


central frequency Iie = OUU0 riz 

FM-modulation — linear excursion of 500 Hz centered at ih 
Sampling frequency eile: = 1000 Hz 

number of samples —N = 256 

pulse duration or integration time —T = N/fs = 256 msec 


Since W was unknown the number of samples (N) could be greater, equal 
or smaller than the number of independent samples (n). However, consid- 
ering that the signal power spectrum is entirely contained in the noise 
spectrum and assuming the noise-spectral-density constant over the physical 
bandwidth B, W = B for the linear case and is greater for the nonlinear 


case. Thus, 


n= 2BT = 2x500x,.256 = 256 
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This was then regarded as the linear processor having all the 
independent samples it was able to process, while the other two had 
less information than they were capable of processing. 

Analysis of the results showed that the linear cross-correlator is 
always better than the nonlinear processors but it is necessary to 
consider that the last two were processing less information than they 
could for the preset bandwidth. 

Also, it was observed that for small signal-to-noise ratio the 
performance of the nonlinear processors approached the linear processor. 
This is explained by inspection of Shannon's expression. As SNR increases, 
the nonlinear processors are unable to process the amount of information 
available since the amplitude limiting allows only two identifiable 
States. This process imposes an upper bound of 2BT bits of information 


as seen by the processor, In other terms, since 


SNR, = f(p,w)SNR, then 
ASNR, = Aftpiu) Aw SNR, + £(p,w)ASNR, 


For the nonlinear processors, when W increases f(p,w) approaches 


f(p) = const and ASNR, = ASNR, x const 


while with the linear correlator f(p,w) increases with W. 

At the lower end, if N is such that the’full information may be 
processed by the nonlinear correlators, the wider bandwidth enhances 
their performance and no gain is achieved by the linear processor, since 


the extra samples are redundant or beyond its processing capability. 


C. INTEGRATION TIME X BANDWIDTH 
In this test the same FM signal was used with a reduced integration 
time. 
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To decrease the pulse duration, the number of samples was reduced 
to N = 128 while the sampling frequency was unchanged. Then T became 
128 ms and the information content was cut to half of the previous value. 

Only the clipper-limited processor was used. Comparison of these 
results with the others previously obtained revealed that the performance, 
or processor, gain as well as the detection sensitivity was rediced as 
expected. 

A similar result could be achieved if the pulse duration was kept 
unchanged and the bandwidth divided by two since in this case the 


information content would also be halved. 


D. INFORMATION BANDWIDTH 

In order to explore further the meaning of the information bandwidth, 
a CW burst code at the Frequency of 5000 Hz was generated. The number 
of samples was again 256 with a sampling frequency of 10 kHz with a 
resultant burst of 25.6 msec. The bandwidth of this pulse was approxi- 
mately 40 Hz. The input filter was kept adjusted for a bandpass of 400 
eee 

A second run was also made with a sampling frequency of 100 kHz, 
reducing T to 2.56 msec and increasing B to 400 Hz. 

The conclusion derived from the analysis shows taht in the first case 
the processor "gain" was much higher than in the second. However, the 
ROC curves show that this contributed nothing to the detection sensitivity, 
but, on the contrary, decreased performance as is shown by the increase of 
false alarm rate. The reason for this rests on the fact that the bandwidth 
used was much higher than the information bandpass, and therefore unneces- 
sary noise power was admitted over a band where the signal presented no 


Significant energy. 
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This confirms the concept that processor gain is a significant 
parameter only if the bandwidth utilized corresponds to that compatible 


to the processor, and to signal and noise characteristics. 
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VI, CONCLUSION AND RECOMMENDATIONS 


The final result of the study revealed the system as a powerful 
instrument for signal analysis. It can be used as an aid to design work 
or as a complement for theoretical studies in signal processing and 
detection. 

The most important feature achieved, out of those initially specified, 
was that of the interaction between man and computer. Although still 
limited to some procedures, the operator can develop freely his reasoning 
in a successive stream of requests. Signal analysis can then be performed 
in a conversational and profitable manner. 

In order to have this capability, the flexibility of the general 
purpose computers was efficiently explored. And no limits exist in this 
respect. Many more situations can be idealized and simulated, and other 
applications devised. 

One field that immediately deserves closer attention is that of actual 
Signal analysis. The system can handle this task without difficulty; 
however, only a continued experience can show which improvements are to 
be made. Depending on the nature of the problem it may be advisable to 
make some changes in order to trade flexibility and versatility for speed 
and precision. In this case, for instance, simplified versions of the 
Subroutines FASTCOR, ANALYSIS and PROCESS could be used without overlay 
and with higher sampling rates to obtain very good statistical averages 
in a reasonable time. 

Some other applications the writer had in mind to verify time permitting, 
are commented upon next, They are a suggestion for future studies with 


the system. 
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Input filter: More elaborate filters either analog or digital may be 
Simulated to represent a transducer accurately. 

Transmission medium: The effect of the transmission medium on the wave 
propagation may be well represented when the characteristics are known. 
For instance, the frequency selective transmission loss may be repre- 
sented by a 20 db per decade type filter and the temperature and 
Salinity influence on the sound velocity could be implemented based 
on the equations giving c = We [Ref. 9]. 

Multipath effects can be analyzed by using as a signal a combination of 
Samples generated according to a predetermined delay pattern. 

Doppier: A very simple way to simulate doppler shift situations is to 
detune the pattern generator during the replica sampling periods. 

Statistical Wave Period Processors: This processor simulation can be 
obtained by a modification in the subroutine FASTCOR to count zero 
crossings of the cross-correlation function and then compute statistics 
upon this data. 

Filter analysis: An example of this possibility was worked with the 
narrow band filter. Other filters either analog or digital may be 
analyzed with no difficulty. If a simple program is added, it would 
be possible to sketch desired filter waveforms on the screen and then 
analyze them. 

Buopeller signature recognition: To treat this very exciting problem 
the harmonic generator may simulate a signal with a well-defined line 
spectrum. It is not however an easy signal to detect when it is 
buried in noise, with low sianaietocneice ratios. The recursive mode 
of the correlation operation is appropriate to perform this computation 
without occupying excessive storage. This should be a good starting 


point to develop a special processor to handle this situation. 
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In conclusion, it is felt that the signal processing system presented 
here can be an effective tool for the exploration of processing methods 


in both instruction and research. 
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APPENDIX A 
SYSTEM OPERATION 


This Appendix includes a description of how to control the system 
from the graphics terminal. 

Before loading the program, it is required to have accomplished 
the items below: 


a) Assign a unit to file 3. 


b) Connect the external clock (synthesizer) to C06 in the analog 
board and start the built-in clock. 


c) Connect the long cable switch between the upper leftmost tie 
point in the logic board and ground. 


d) Load and execute GATED in the AGT. 

After loading, the teletype prints a reminder to have the analog 
ready, AGI GATED and clock running and asks for ID, the terminal 
identification number. When the message ID=l(or 2} is entered, the 
operation starts. 

At the terminal only the keyboard is necessary to type the commands. 

In the beginning, the constants are displayed and can be modified 
by erasing the block and typing the new value with the first digit 
aligned with the first letter of the heading. After carriage return the 
SETPOT mode is entered to initialize potentiometers. Then, the system 
waits for the first command. 

Whenever an illegal operation or value is typed, the message 
“Incorrect Selection, Try Again" is displayed for a few seconds. 

The commands are as follows: 

A. CONSTANTS 
M - Sample block size; initial 256; may be any power of two less 


than or equal to 256 and greater than 8. 
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SAMP 


THRS 


LINE 


PROC 


POIKS 


Sep 


CPDT 


LNTH 


1, 2, or 3 indicates to sample replica, signal initial, or 
Signal continuation; initialized with 1. 

threshold setting; 1000 means 1.0 normalized threshold; greater 
values are illegal; initialized with 0. 

line number of the sampling point; may be any ADK legal calling 
number; initialized with 1. 

type of processor; 0 is linear processor, 4 is clipper limited 
and 8 is amplitude limited; initialized with 0. 

initial value of pattern generator injection potentiometer; 
initialized with 9999; the noise generator potentiometer value 
is the complement to 10000; in a manual operation needs to be 
set only if the required sequence is correlation or convolution 
recursive. 

value of decreasing step for signal injection potentiometer to 
be used in performance or ROC determination; must be less than 
POTS; initialized to 4000. 

if 5 indicates that false target rate is required; if 4 | 
initiates a distance determination; if 3 a probability of 
detection; these operations, however take place only if a 
correlation recursive is asked for next; initialized to one. 
length parameter for LPLOT; has to be an even number between 


4 and 16; initialized to 16 (smaller plot size). 


All these are four-digit numbers and modifiable. The next three are 


displayed only. 


SPOW - 
NPOW - 


input signal power; initial 1.0 


input noise power; initial 0 


SNR (DB) - signal to noise ratio in db, initialized to 0. 
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B. ANALOG CONTROL 

SETPOT - displays alternately the words "POT VAL"; allowable 
Pie eieneten identifiers and values may be typed left justified under 
the first letter; sets the potentiometers. 

SETLINES = displays alternately "SET VAL"; allowable line identifiers 
(two digits) and ] or -1 may be typed left justified; sets the lines. 

TEST - displays alternately "TST VAL"; allowable line identifiers 
may be typed; after return displays | or -] according to the test; 
requires another return. 

POTSET, COMPUTE, HOLD - puts the analog in the mode specified. 

SCAN - displays current values of all amplifiers, potentiometers and 
several other active components, read in the current analog mode; waits 
for the message REIN to return or POTS to scan the potentiometers in 
POTSET mode; if the last option was entered only the potentiometer values 


are displayed and RETN must be typed to return. 


C. ANALYSIS MODE 

SAMP - asks for a confirmation about the constant values; if NO is 
typed, it brings in CONSTANTS; if return is keyed, executes sampling. 
(If a recursive correlation or convolution is to be performed, the 
replica must be sampled first with the potentiometer set to 9000 and next 
the signal, with the signal input potentiometer adjusted to the value 
read in POTS). 

FFTR - asks for function to be transformed and displays function, 
absolute value and phase of the transform; results are not stored. 

CORRELATION - asks for type and displays replica, signal and 


correlation; if recursive, the operation is terminated by the operator 
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closing the manual switch for a moment; also shows the threshold line 
and information about threshold reached, probability of detection 
(meaningful only if CPDT were previously set to 3 or 5), signal-to-noise 
ratio (the last one computed) and distances (if CPDT is set to 4). 

CONVOLUTION - asks function to be convoluted with the replica, and 
type; displays replica flipped, the other function and convolution; if 
recursive, stopped by the manual switch. 

POWER - asks for function and type (recursive not allows); displays 
autocorrelation and power spectrum (absolute value and phase); results 
ave ot stored. 

AUTOCORRELATION - asks for function and displays function (twice) 
and result. 

DISPLAY - asks for function and displays. 

SNRI - external clock must be set for time scale operation; result 


may be obtained through CONSTANTS. 


D. SIMULATION MODE 

CWBURST - adjusts the circuits for a preset CW burst code 

FMSLIDE - adjusts the circuits for a preset linear FM slide code. 

FSK - sam for FSK code. 

PROPELLER - same for a propeller signature; no automatic operations 
were implemented for this signal. 

PERFORMANCE - computes the input-output characteristics for the 
processor and signal selected, and displays the result; external clock 
must be set for time-scale; an extra return required to end operation. 

ROC - same as above but with the computation of Receiver Operating 
Characteristic 


DIST - presets the circuits and initiates distance measurement operation. 
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few PLOT 

VPLOT - plots a VPLOT on the line printer. 

LPLOT - plots a LONGPLOT on the line printer with dimension 
inversely proportional to LNTH. 

RPLOT - asks to ready the recorder and informs the setting for 
the probability counter potentiometer if necessary; after return, informs 
that a delay of 20 seconds is in course to allow starting the eight- 
channel recorder; during this period a zero is sent to channels 6, 7 
and 8 for calibration; the amplitude scaling is to be 2v per division 


and the speed 5 or 10 mm. per second. 


Gee OP 


Ends the program and releases the digital computer. 
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APPENDIX B 
ANALOG COMPUTER INFORMATION 


In this Appendix are derived the logic equations to control the 
PATTERN GENERATOR and the FSK SYNCHRONIZER, Also are included tables 
for quick scaling the potentiometers, (TABLE I), for determination of 
the coefficients of the frequencies generated by the HARMONIC GENERATOR 
(TABLE II) and for adjustment of the shift frequency of the FSR and the 
low pass RC filter (TABLE III). 

The following variables are defined: 


RST - general reset mode 

CMP - general compute mode 

HLD - general hold mode 

RRP - Ramp reset mode 

CRP - ramp compute mode 

HRP - ramp hold mode 

RO6 - switch ROO NC 

R10 - switch R10 NC 

R26 - switch R26 NC 

R25 - R input of A025 

H25 - H input of A025 

T07 - line T007 (program call 8) 
T10 - line TO10 (program call 9) 
T11 - line TO11 (program call 10 
T15 - line TO16 (program call 15) 
T17 - line T017 (program call 16) 
C04 - comparator C04 

CO7 - comparator CO/ 

CWB - CW burst 

FM - FM slide burst 

FSK - FSK burst 

PRO - propeller simulation 
DIST - distance burst 


Ps 


PTS - potset 

SHC - clock line (FSR-shift) 
RLN - reset line (FSR) 

CDV - frequency divider 

SIN - sample interrupt enable 


PATTERN GENERATOR 
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RO6 = TI7 = DIST + PRO 
R16 = 116 = FSK 
RST = T07 - T10 
CMP = TO7 - TI10 
HLD = TO7 - T10 


The above equations are controlled by the program exclusively through 


CALL SETLINES statements. 


feo = HRP + CRP 


H25 = RRP + CRP 


| 
| 
| 
| 
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HRP = HLD - CO7 - CWB = CW CO7 - TO7 - TI10 

CRP = CMP - CO7 : CWB = CWB - C07 - TO7 - T10 

RRP = RST + CO7 + CWB = CWB + CO7 + TO7 - T10 

R25 = CWB - CO7 - TIO - TO7 + CWB - CO7 - T10 - TO7 
erp CO, - 110 = CHR - G07 - Td 
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and then 


R25 = CWB 
= CWB 


Making CW 


R25 = Tl] - 


= Tl] also controlled by the program, the equations: 


ed 


O7 - T10 3; H25 = T11 + TO7 and RIO = C04 were implemented 


in the logic board as shown in Fig. 3.3. To prevent possible occurrence 


of a race 


problem and oscillations a short duration latch was included. 


SYNCHRONIZER 


SHC = 


SIN = 
RLN 


a eee ee 


CMP + PTS + CDV = TO7 - T10 + PTS + CDV 


| 


a i 
EN a a 


Oye 110 PTS + COV ee O7 cee eee See 


© 
© 
<= 


ee 
re —=-— —_—- ———— —_——_— 


= 107 - [Oe eas 


“5 
YN 


Chas 


S| 


—— 


= T07 - T10 : POT = TO7 +: T10 : PTS 


| 
| 


0 
Sal 


0 


The implementation is shown in Fig. 3.9. 
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TABLE I 
POTENTIOMETER SETTING 


{| He 


Real time C= 19 
Time scale C= .0l 
Ea oe fsc = ‘ PS fe. 

F P F P E P E P 
1.59 9999 1.19 7480 79 4960 39 2445 
ee O7 9986 hay, F350 1] 4840 mi) Does 
a0 8874 els 7230 5 4715 35 2198 
le 53 8861 | Bs 7200 ay es 4585 33 2055 
ei) 9495 ea 6970 el 4455 230 1948 
1.49 9352 1.09 6840 .69 4340 29 1821 
1.47 C285 07 6720 6/7 4205 27 1695 
1.45 9120 | 50s 6595 265 4085 OE 1570 
1.43 8995 0s 6465 51a) 3948 123 1443 
1.41 8860 1.01 6345 61 3835 21 1319 
ire? S725 .99 6220 38) 3703 19 1193 
Meo 7 8610 97 6095 Sy 3580 ay 1058 
eS 8490 5 5965 oo 3478 15 0942 
33 8351 93 5840 og aa 13 0816 
te 3) 8230 91 5720 “e)| 3205 mi 0691 
lee? 8102 89 5590 49 3078 .09 0565 
t.2/ 7970 87) 5460 ay 9935 var 0439 
e225 jisseyl 85 5340 »45 2095 05 0314 
23 7730 ua: aye |l8 43 2700 n0G 0188 
Wee | 7600 .8] 5095 4] 2575 50) 0063 





TABBES EI 
HARMONIC CONTENT 


e= f(f)) 
Frequency Coefficient 
2 4 poe 
DC = Kna /2 + 3K,a /8 - Skea /16 
3 5 
f Kya - 3K,.a Vil as oKea /8 
2 4 6 
of Kna /2 - K,a /2 + 15Kea V2 
3 5 
she he - Kea /16 
6 
Af K,a /8 - 3Kea /16 
6f Kea” /32 
e = F(f,, fy) 
Order Frequency Amp. due to lowest K 
] Ky oK skp f, »f, Kya, K,b 
2 
2 Kn Ka ke 2f, »2f, 1/2 Koa 
f, +f, Koa 
3 
3 Ko Ke 3f, »3f, 1/4 Ee 
2f,+f, e/a Ka b 
4 
4 Ky Ke Af, 4, 1/8 Ky? 
3f, +f, 1/2 Ka b 
2f +f. 
5 K Sf, ,5¢ 1/16 Kea” 
5 l2=22 5 A 
Af, +f, 5/16 Ka b 
3f,#2f, 
6 K 6, ,6f 1/32 Ka 
6 | aa 6 5 
of) +f. 3/16 Kea b 
Af, +2f, 
3f,#3f, 


117 





TABLE III 
FSR FREQUENCY AND RC FILTER 
f 





fcs = oS 
Built in Clock Divider + _ fcs R C Series/Paral lel 

100 8 G25 Al 2x10 S 
3x] p 
100 4 2.0 “1 2x10 S 
100 2 3.38 aI 3x] p 
100 5 1 2x] Pp 
1000 8 1255 al 2x1 S 
oad P 
1000 4 20 © (| 2x] S 
1000 2 3355 d Sy] P 
1000 1 50 sl Axa 2 
10000 8 125 3x] S 
10000 4 200 al 2x10 S 
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COMPUTER PROGRAM 
FOR 
HYBRID SONAR SIMULATION 


MAIN PROGRAM 

Purpose: The main program is used for computer initialization 
and serves as a bridge to keep the status during overlay swapping. It 
also creates two dummy pulses as Replica and Signal for operative 
checking after initialization, | 

Parameters and variables: In the explanation of the program the 
meaning of a variable is given only in the subroutine where it is 
effectively used. In the case of the MAIN program the variables below 
are of interest: 


M- basic sample size. Has to be a power of two. Initially has the 
value 256. 


YR - array for storage of Replica, within size M. Initially 
dimensioned to 256. 


YS - array for storage of signal, with size 3M. Initially dimensioned 
to 768. 


YC - array for storage of correlation, with size 2M. Dimensioned 
Om otc. 


VAGO - buffer, with size 5M and dimensioned to 1280. 
IB - buffer, with size 2M plus one word. Dimensioned to 1025. 
MI - text directory. 
MG - graph directory. 


ID - graphics terminal number. 
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